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Abstract. Forty years elapsed since the optical identification of the first isolated neutron star (INS), the Crab pulsar. 25 
INSs have been now identified in the optical (O), near-ultraviolet (nUV), or near-infrared (nIR), hereafter UVOIR, including 
rotation-powered pulsars (RPPs), magnetars, and X-ray-dim INSs (XDINSs), while deep investigations have been carried out 
for compact central objects (CCDs), Rotating RAdio transients (RRATs), and high-magnetic field radio pulsars (HBRPs). In 
this review I describe the status of UVOIR observations of INSs, their emission properties, and I present recent results. 
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INTRODUCTION 

Interestingly enough, although isolated neutron stars 
(INSs) were discovered in the radio band as sources of 
strongly beamed radiation (Hewish et al. 1968), hence 
dubbed pulsars, one object of this class had been already 
observed for many years back in the optical band. This is 
the Baade's star, a quite bright object {V ~ 16.6) located 
at the centre of the Crab Nebula in the Taurus constella- 
tion. It was only after the discovery of a bright radio pul- 
sar (NP 0532) at the centre of the Crab Nebula (Comella 
et al. 1969) that the Baade's star was indeed recognised 
to be a possible INS and certified as such by the discov- 
ery of optical pulsations at the radio period (Cocke et al. 
1969). 

Almost a decade passed before the optical identifica- 
tion of another INS: the Vela pulsar (PSR B0833-45), 
one of the faintest objects ever detected on photographic 
plates (V = 23.6; Lasker 1976), soon after identified 
as an optical pulsar (Wallace et al. 1977). While the 
search for optical counterparts of radio pulsars continued 
through the years and scored several new identifications 
(see, e.g. Mignani et al. 2004 for a review), other classes 
of typically radio-silent INSs were progressively discov- 
ered at X-ray and y-ray energies, which also become in- 
teresting targets for optical, as well as near-ultraviolet 
(UV) and near-infrared (IR) observations. These are the 
Soft Gamma-ray Repeaters (SGRs) andf the Anoma- 
lous X-ray Pulsars (AXPs), the magnetar candidates (see 
Mereghetti 2008; 2009), the X-ray Dim INSs (XDINSs; 
Haberl 2007; Triimper 2009), and the central compact 
objects (CCOs) in SNRs (De Luca 2008). At variance 
with radio pulsars which are powered by the neutron star 
rotation (Gold 1968; Pacini 1968), and then are usually 
referred as rotation-powered pulsars (RPPs), these INSs 
radiates through different emission mechanisms like, e.g. 
the decay of an hyper-strong magnetic field, for the mag- 
netars, or the cooling of the neutron star surface, for 



the XDINSs. At the same time, peculiar radio-loud INSs 
were also discovered, like the high-magnetic field radio 
pulsars (HBRPs; Camilo et al. 2000) and the Rotating 
RAdio Transients (RRATs; Mc Laughlin 2009). In many 
cases, optical or nIR counterparts have been identified 
both for the XDINSs and for the magnetars. 

Thus, despite of their intrinsic faintness which makes 
them very elusive targets, nUV, optical, and nIR (here- 
after UVOIR) observations proved quite successful in 
detecting INSs outside the radio band, thanks to both 
the HST and to the new 8m-class telescopes, like the 
VLT. Forty years after the optical identification of the 
Crab pulsar, UVOIR studies of INSs represent an impor- 
tant tile to complete an understand the multi-wavelength 
phenomenology of these objects. In this paper, I will re- 
view the status of the UVOIR observations of all differ- 
ent classes of INSs, their emission properties, and the 
results of recent observations. 



RPPS 

Being the first class of INSs detected at UVOIR wave- 
lengths, the major observational effort has been devoted 
so far to the search and to the study of RPPs which, in- 
deed, score the largest number of INS identifications. 

After the identification of the Crab and Vela pulsars, 
the major advances in the UVOIR studies of INSs were 
achieved thanks to observations performed with the ESO 
NTT 'm the late 1980s and in the 1990s (see, e.g. Mignani 
et al. 2000 for a review). These yielded to the iden- 
tification of the optical counterparts to the first radio- 
silent INS, Geminga (Bignami et al. 1993), to the first 
extra-galactic INS, PSR B0540-69 in the LMC (Car- 
aveo et al. 1992), and to PSR B0656H-14 (Caraveo et 
al. 1994). The refurbishment of the HST in 1993 rep- 
resented a major turn-off in UVOIR studies of INSs 
(see, e.g. Mignani 2007), yielding to the nUV identifi- 



cation of PSR B0950+08, PSR B1929+10 (Pavlov et 
al. 1996), PSR B1055-52 (Mignani et al. 1997), and 
of PSR J0437-4715, the first ms-pulsar detected at 
UVOIR wavelengths (Kargaltzev et al. 2004). Moreover, 
HST observations allowed to study, for the first time, 
the morphology and evolution of pulsar-wind nebulae 
(PWNe) around the Crab pulsar (Hester et al. 1995) and 
PSR B0540-69 (Caraveo et al. 2000). In the 2000s, ob- 
servations with the VLT yielded to the identification of 
the optical counterparts to PSR B1509-58 (Wagner& 
Seifert 2000) and, more recently, to PSR Bl 133+16 
(Zharikovetal. 2008) and PSR JO 108 -1431 (Mignani et 
al. 2008a). A possible counterpart to PSR J0108-1431 
was indeed spotted a few years before by Mignani et al. 
(2003), barely detected against the halo of a nearby ellip- 
tical galaxy, but it was not recognised as such until recent 
Chandra observations discovered X-ray emission from 
the pulsar (Pavlov et al. 2009) and measured its proper 
motion {p. = 200 ±65 mas yr^') through the compari- 
son with the radio position. This allowed Mignani et al. 
(2008a) to find that the position of the candidate counter- 
part (U = 26.4) nicely fitted the backward proper motion 
extrapolation of the pulsar and, thus, to certify its identi- 
fication a posteriori 

For many objects, the faintness of the optical counter- 
part (V > 25) initially prevented the use of the straight- 
forward optical timing identification technique, success- 
fully used for the brighter Crab and Vela pulsars and for 
PSR B0540-69 (Shearer et al. 1994). In some cases, 
the proper motion measurement of the candidate coun- 
terpart, successfully tested for Geminga (Bignami et al. 
1993), was used instead, confirming the identification 
of PSR B0656+14 (Mignani et al. 2000) and of PSR 
B 1929+10 (Mignani et al. 2002), both achieved through 
high-resolution HST imaging. For both Geminga and 
PSR B0656+14, the identification was later strength- 
ened by the detection of nUV pulsations with the HST 
(Kargaltsev & Pavlov 2007). For the others, the identi- 
fication evidence mainly reUes on the positional coin- 
cidence with the INS coordinates (PSR Bl 133+16 and 
PSR J0108— 1431), on the peculiar colours or spec- 
trum of the candidate counterpart (PSR B0950+08), 
on the possible evidence of optical polarisation (PSR 
B1509-58). 

Recently, the optical identification of PSR B 1055— 52 
was confirmed through new HST observations (Mignani 
et al. 2009a). The candidate counterpart was clearly de- 
tected in the nUV with the ACS/SBC (Fig. 1) and in the 
optical with the WFPC2. Together with the original U- 
band HST/FOC photometry of Mignani et al. (1997), the 
new HST flux measurements allowed to measure the ob- 
ject peculiar colours which virtually certify it as the ac- 
tual counterpart of PSR B 1055 —52. Moreover, the com- 
parison between the relative position of the counterpart, 
as measured in the 1996 FOC observations and in the 
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FIGURE 1. HST ACS/SBC image of the PSR B1055-52 
field (12" X 12"). North to the top. East to the left. The pulsar 
counterpart is marked by the arrow. The star ~ 4" north of the 
pulsar is star A of Mignani et al. (1997). 



new ACS ones, allowed to measure a proper motion of 
~ 55 mas yr~\ the first ever measured for this pulsar at 
any wavelength. This corresponds to a tangential velocity 
vt ~ ISOdjQQ, where djoo is the source distance in units 
of 700 pc. Thus, out of the - 1800 RPPs listed in the 
updated ATNF radio catalogue, twelve have been identi- 
fied at UVOIR wavelengths, albeit with different levels 
of confidence. 

The study of the UVOIR emission properties of RPPs 
is complicated by the paucity of spectral information. 
Only four RPPs have a complete UVOIR spectral cov- 
erage, with the nUV and nlR spectral bands crucial to 
disentangle the contributions of thermal and non-thermal 
emission processes. Moreover, only for five RPPs low- 
resolution spectroscopy is available, while for the others 
the knowledge of the spectral energy distribution (SED) 
still reUes on, sometimes sparse, multi-band photome- 
try, with all the implied caveats (see, e.g. discussion in 
Mignani et al. 2007a). Only five RPPs pulsate in the 
optical/UV, where they all feature double-peaked hght 
curves, with the exception of PSR B0540— 69, and ro- 
bust optical polarisation measurements exist only for the 
Crab pulsar (see, e.g. Slowikowska et al. 2009 and refer- 
ences therein), reducing the impact of these important di- 
agnostic tools. The SED of the RPPs grows in complex- 
ity with the spin-down age (see, e.g. Fig. 2 of Mignani 
et al. 2007a). It evolves from a single power-law (PL) 
spectrum (Fy oc v~"; a ~ 0) for young RPPs (age< 10^ 
years) to a composite one featuring both a PL (a > 0) 
and blackbody (BB) component (r - 2 - 8 x 10^ K) for 
the middle-aged ones (age 10^ — 10^ years), the former 



dominating in the IR, the latter in the UV. A similar com- 
posite SED was measured also for the middle-aged PSR 
B 1055 -52 (Mignani et al. 2009a). Surprisingly, older 
RPPs (age> 10^ years) seem to feature single PL spectra. 
For the very old (lO'^ years) PSR J0437-4715, the nUV 
spectrum is consistent with a pure BB (Kargaltsev et al. 
2004). The PL component is ascribed to non-thermal ra- 
diation produced by relativistic particles accelerated in 
the neutron star magnetosphere (e.g., Pacini & Salvati 
1983). On the other hand, the BB component is ascribed 
to thermal emission produced from the cooling neutron 
star surface. In general, there is a correlation between the 
optical luminosity and the neutron star rotational energy 
loss (Zharikov et al. 2006), suggesting that the magne- 
tospheric emission is powered by the star rotation and 
contributes most to the optical luminosity. Emission ef- 
ficiencies Lgpt/E « 10^^ — 10^^ are then derived. It is 
interesting to note that the BB/PL components which fit 
the RPP optical spectra are not consistent with the extrap- 
olation of the analogue BB/PL components which fit the 
X-ray spectra (Mignani et al. 2004). The presence of two 
PL components clearly indicates a break in the magne- 
tospheric spectrum. No evident correlation is found be- 
tween the optical PL spectral index and the neutron star 
age (Mignani et al. 2007a), like none is found in the X- 
rays (Becker et al. 2009). On the other hand, the pres- 
ence of distinct optical and X-ray BB components sug- 
gests that the temperature distribution on the neutron star 
surface is not homogeneous. 

Interestingly, Kargaltsev et al. (2004) noted that the 
BB temperature (J ~ 2 x 10^ K) derived for PSR 
J0437— 4715 is larger than expected from standard cool- 
ing models. This suggests that either the neutron star 
surface temperature was raised through the effect of re- 
heating processes occurring in its interior (e.g., Page 
2009; Tsuruta 2009) or it was raised by accretion from 
its companion star during a past pulsar spin-up phase. 
Like for the old ms-pulsar PSR J0437-4715, re-heating 
might have occurred also for the 200 million year old 
PSR J0108— 1431. The accurate measurement of its dis- 
tance (240 pc), obtained through radio parallax mea- 
surements (Deller et al. 2009), implies that the optical 
fluxes could be compatible with thermal emission from 
the bulk of the neutron star surface at a BB temperature 
of f» 3 X 10^ K. Again, this temperature is larger than the 
values expected from standard cooling models. 

XDINSS 

Soon after their discovery in the 1990s, XDINSs were 
the next class of INSs to attract interest from the neutron 
star optical astronomy community. Indeed, without the 
evidence coming from the detection of X-ray pulsations, 
detected only later with XMM-Newton (e.g., Haberl et al. 



2007) , the measurement of an extreme X-ray-to-optical 
flux ratios Fx/Fopt was, at that time, the only way to 
certify these sources as INSs. Although no XDINS op- 
tical counterpart was initially identified from NTT ob- 
servations, despite of a considerable observational ef- 
fort, they were nonetheless crucial to pave the way to 
deeper observations. The first XDINS to be detected in 
the optical was RX J1856.5— 3754 which was indeed ob- 
served by the //5r (Walter & Matthews 1997). Riding the 
wave, HST also detected likely optical counterparts for 
RX J1308.6H-2127andRX J 1 605. 3H-3249 (Kaplan et al. 
2002, 2003). From the ground, optical detections were 
obtained for RX J0720.4-3 125 with the NTT (Motch & 
Haberl 1998) and with the Keck (Kulkarni & van Kerk- 
wijk 1998) and for RBS 1774 wifli flie VLT (Zane et al. 

2008) and the LBT (Schwope et al. 2009). However, the 
optical identifications were certified by proper motion 
measurements only for three of them: RX J1856.5— 3754 
(Walter et al. 2001), RX J0720.4-3125 (Motch et al. 
2003), and RX J1605.3H-3249 (Motch et al. 2005; Zane 
et al. 2006), while for botii RX J1308.6h-2127 and RBS 
1774 the optical identifications are still based on the 
positional coincidence with the Chandra X-ray coordi- 
nates and on the Fx/Fopt ratio. For the two XDINSs for 
which optical parallaxes were measured with the HST, 
RX J1856.5-3754 and RX J0720.4-3125 (van Kerk- 
wijk & Kaplan 2007), the proper motion measurement 
also allowed to derive their space velocities and, thus, to 
rule out accretion from the interstellar medium (ISM). 

Recently, a possible candidate counterpart was identi- 
fied for RX J0420.0-5022. OriginaUy, a candidate coun- 
terpart was tentatively proposed by Haberl et al. (2004), 
based on VLT observations. However, a re-analyses of 
their VLT data (Mignani et al. 2009b) showed fliat flieir 
candidate was detected only at the ~ 2a level and it was 
most likely a spurious detection produced by the halo of 
a very bright {B= 10) star located 40" away. At the same 
time, deeper VLT observations (Mignani et al. 2009b) 
allowed to find evidence for the detection (~ 4(7) of a 
fainter object {B = 27.5 ± 0.3) coincident with the Chan- 
dra position of RX J0420.0-5022 (Fig. 2). 

XDINSs are only detected at optical wavelengths. 
Although they have been repeatedly observed in the 
nIR with the VLT, none of them was detected so far 
(Mignani et al. 2007b, 2008b; Lo Curto et al. 2007; 
Posselt et al. 2009). No observation of XDINSs in tiie 
nUV was reported yet. Like for the RPPs, the knowl- 
edge of the SED mainly rehes on multi-band photom- 
etry, with low-resolution spectroscopy only performed 
for RX J1856.5-3754 (van Kerkwijk & Kulkarni 2001). 
No search for optical pulsations or polarisation measure- 
ments have been carried out so far. 

The XDINS optical fluxes usually exceed by a factor 
of ~ 5 (or more) the extrapolation of the X-ray black- 
body (see, e.g. Kaplan 2008). For the XDINSs wifli an in- 




FIGURE 2. Left: VLTIF0RS2 B-band image (1' x 1') of the RX J0420.0-5022 field. North to the top, east to the left. The circle 
(1.11" radius; 90 % confidence level) indicates the RX J0420.0-5022 position (Mignani et al. 2009b). A bright star south of RX 
J0420.0— 5022 is masked. Right: 10" x 10" zoom of the same region smoothed with a Gaussian filter over cells of 5 x 5 pixels. 
Object X (S = 27.5 ± 0.3) is the putative RX J0420.0-5022 counterpart. Field objects are labelled as in Haberl et al. (2004). 



ferred rotational energy loss, i.e. RX J0720.4-3125, RX 
J1308.6+2127, RX J1856.5-3754, RBS 1774 (Kaplan 
& van Kerkwijk 2005a,b; van Kerkwijk & Kaplan 2008; 
Kaplan & van Kerkwijk 2009) this turns out to be too 
low (« 10^" erg s^') to sustain magnetospheric emission 
powered by the star rotation, unless the emission effi- 
ciency is a few orders of magnitude higher than in RPPs. 
In at least in the best-studied cases of RX J1856.5— 3754 
(van Kerkwijk & Kulkarni 2001) and RX J0720.4-3 125 
(Motch et al. 2003), the optical spectra is found to closely 
follow a Rayleigh-Jeans, which suggests that the optical 
emission is predominantly thermal (T ~ 2 — 4 x 10^ K). 
The XDINS optical emission has been thus interpreted 
either in terms of a non homogeneous surface tempera- 
ture distribution, with the larger and cooler part emitting 
the optical (e.g.. Pons et al. 2002), or of reprocessing of 
the surface radiation by a thin H atmosphere around a 
bare neutron star (Zane et al. 2004; Ho 2007). For RX 
J1605.3H-3249, multi-band photometry measurements do 
exist (Motch et al. 2005) but they are likely affected 
by instrument cross-calibration problems which hamper 
the SED characterisation (Zane et al. 2006). For RXJ 
1308.6H-2127 (Kaplan et al. 2002), RBS 1774 (Zane et 
al. 2008; Schwope et al. 2009), and RX J0420.0-5022 
(Mignani et al. 2009b) only one band measurement is 
available 

Interestingly, RBS 1774 is characterised by an anoma- 
lously large optical excess of « 35 (Zane et al. 2008). 
This is unlikely produced from the cooling neutron star 
surface unless, to explain the low observed pulsed frac- 
tion of the X-ray light curve, one invokes a peculiar 



alignment between the magnetic field, the rotations axis, 
and the fine of sight (Zane et al. 2008). The RBS 1774 
emission may thus be of magnetospheric origin, perhaps 
related to its large magnetic field (~ lO'"*^ G), inferred 
from the observations of an absorption feature in the X- 
ray spectrum (Zane et al. 2005). However, the recent 
measurement of the pulsar spin-down (Kaplan & van 
Kerkwijk 2009) implies a magnetic field of 2 x lO'^ 
G, thus possibly arguing against this interpretation. A re- 
cent re-analysis of the XMM-Newton spectrum of RBS 
1774 using a different spectral model (Schwope et al. 
2009), while confirming the large optical excess, does 
not rule out the possibility that the optical emission is 
indeed thermal. For RX J0420.0-5022, the flux of the 
putative counterpart coiTesponds to an optical excess of 
^ 7 with respect to the extrapolation of the XMM-Newton 
spectrum (Mignani et al. 2009b). From the upper limit 
on the RX J0420.0-5022 rotational energy loss, E < 
8.8 X lO-'^ erg s^' derived from X-ray timing (Haberl et 
al. 2007), the optical luminosity of the putative counter- 
part. Lb ^ 1.2 X 10^^ erg s^' d^^Q (where d$5o is the dis- 
tance in units of 350 pc; Posselt et al. 2007), implies an 
emission efficiency > 1.3 x 10^^ for rotation-powered 
emission. This value could still be compatible with the 
range of emission efficiencies derived for 10*' — 10^ years 
old RPPs (Zharikov et al. 2006). In case of thermal emis- 
sion from the neutron star surface, the optical luminosity 
of the putative counterpart would be compatible with a 
blackbody with temperature < 25 eV and an an implau- 
sibly large emitting radius of > 23 km. Thus, would the 
putative counterpart be confirmed, its optical luminosity. 



for a neutron star distance of ^ 350 pc, might rather point 
towards a non-thermal origin for the optical emission. 

SGRS AND AXPS 

Optical/IR observations of AXPs and SGRs, the magne- 
tar candidates, were originally carried out (e.g., Davies 
& Coe 1991; van Paradijs et al. 1996) to investigate the 
nature of these sources and to test models based on ac- 
cretion from a companion star or a debris disc. Unfortu- 
nately, optical/nUV observations of magnetars are ham- 
pered by the high interstellar extinction (up to Ay ~ 30) 
in the direction of the galactic plane, where most of them 
were discovered. This makes the nIR the only spectral 
band suitable for observations, although they are compli- 
cated by the field crowding. Indeed, nIR identifications 
of magnetars were boosted when deep adaptive optics 
(AO) observations became possible with 8m-class tele- 
scopes like the VLT, Gemini, as well as with the CFHT, 
allowing to perform high-resolution imaging and to re- 
solve potential candidate counterparts within the Chan- 
dra X-ray source error circles. However, unless a more 
accurate source position is available from radio obser- 
vations like, e.g. in the case of the two transient radio 
AXPs XTE J1810-197 and IE 1540.0-5408 (Camilo 
et al. 2007a,b) and of SGR 1806-20 (Cameron et al. 
2005), the positional coincidence alone is not sufficient 
to pinpoint the actual counterpart. Given the source dis- 
tance and the longer required time span, proper motion 
measurements are not a very efficient identification tool 
Uke, e.g. in the case of RPPs and XDINSs, while the 
search for nlR pulsations at the X-ray period requires 
suitable instruments which are usually not available at 
most telescopes. Indeed, optical pulsations have been de- 
tected from two AXPs (see below) using guest instru- 
ment facilities. Moreover, the still significant interstellar 
extinction, together with the uncertainty on the source 
distance and the difficulty of finding a well-defined spec- 
tral template for magnetars, makes a colour-based identi- 
fication uncertain. Thus, in most cases the best recipe for 
nIR identifications of magnetars is based on the measure- 
ment of correlated IR/X-ray variability. The application 
of this recipe, of course, benefited from the possibility 
of performing prompt Target of Opportunity (ToO) ob- 
servations with ground-based telescopes in response of 
triggers from high-energy satellites. 

Out of the five currently known SGRs, SGR 1806-20 
has been the first to be firmly identified in the nIR (Is- 
rael et al. 2005; Kosugi et al. 2005). Recently, a sec- 
ond nlR identification has been obtained for the newly 
discovered SGR 0501H-4516 thanks to prompt follow-up 
observations with the UKIRT, the WHT, and the Gemini 
(Rol et al. 2009) which pinpointed a ^Ty ~ 19.1 coun- 
terpart at the boresight-corrected 0.1" Chandra position 



(Woods et al. 2008). The same observations also showed 
a possible evidence for nlR variability, although not evi- 
dently correlated with that observed in the X-ray band. 
The counterpart was also detected in the optical (see 
Rol et al. 2009 et references therein). At the same time, 
four out of the eight known AXPs have been identi- 
fied in the nIR: 4U 0142H-61 (HuUeman et al. 2004), 
IE 1048.1-5937 (Wang & Chakrabarty 2002; Israel et 
al. 2002), XTE J1810-197 (Israel et al. 2004; Rea et 
al. 2004), and IE 2259H-586 (Hulleman et al. 2001). 
In addition, both 4U 0142H-61 (Hulleman et al. 2000) 
and IE 1048.1-5937 (Durant & van Kerkwijk 2005) 
have been identified in the optical. 4U 0142H-61 and 
IE 2259H-586 are the only AXPs/SGRs which has been 
detected in the mid-IR by Spitzer (Wang et al. 2006; 
Kaplan et al. 2009), while upper limits have been ob- 
tained for XTE J1810-197, IRXS J170849.0-400910 
(Wang et al. 2007a) and IE 1048.1-5937 (Wang et al. 
2007a; 2008). nIR counterparts were also proposed for 
SGR 1900H-14 and for the AXP IE 1841-045 (Testa 
et al. 2008), based on possible long-term nIR variabil- 
ity. Optical flares of a possible candidate SGR, SWIFT 
J195509.6H-261406, were detected by Stefanescu et al. 
(2008) using the OPTIMA instrument at the 1.3-m tele- 
scope of the Skinakas Observatory (Crete), with the pos- 
sible evidence of periodicity (6-8s) in the strongest flares. 

Recently, an nlR counterpart was identified for the 
radio-transient AXP IE 1540.0-5408 flirough VLT ob- 
servations performed in January 2009, right after the 
source underwent an X-ray outburst detected by the 
Swift/BAT (January 22; GronwaU et al. 2009). A candi- 
date counterpart (Fig. 3) was identified coincident with 
the radio error circle of IE 1540.0—5408 (Camilo et al. 
2007b) with a magnitude = 18.5 (Israel et al. 2009a). 
The candidate counterpart was not detected in Magellan 
observations (Ks > 20.1) obtained just four days before 
the X-ray outburst (Wang et al. 2009). The candidate 
counterpart was also not detected in much deeper VLT 
images of the same field acquired in July 2007 (Mignani 
et al. 2009c) down to a Umiting magnitude of Ks ~ 21.7. 
This implies a long-term variability of about 3 magni- 
tudes (Kg band), which virtually certifies its identification 
wifli IE 1540.0-5408. 

It is interesting to note how all magnetar identifi- 
cations were achieved with ground-based telescopes, 
mostly with the VLT, while the HST did not play a sig- 
nificant role, partially because its nIR camera (NICMOS) 
was unavailable between 1999 and 2002. 

Historically, optical/nIR observations of SGRs and 
AXPs were fundamental to provide crucial supporting 
evidence for the INS scenario and for the magnetar 
model. In the IR, magnetars are characterised by a flux 
excess with respect to the extrapolation of the X-ray 
spectrum (see, e.g. Israel et al. 2005). Although this has 
been often referred in the Uterature to be a distinctive 



FIGURE 3. VLTINACO band images (10" x 10") of the AXP IE 1540.0-5408 field obtained in July 2007 (left; Mignani et al. 
2009c) and in January 2009 (right; Israel et al. 2009a). North to top, East to the left. The variable object in the computed 3<T radio 
error circle (0.63") is the nIR counterpart. 



character of magnetars, many RPPs indeed feature a sim- 
ilar excess (see, e.g. Mignani et al. 2007a) and it can be 
interpreted in terms of either a spectral break or of the 
onset of a new spectral component. 

The magnetar optical/nIR SED is not fully charac- 
terised yet. Indeed, in most cases the accessible spec- 
tral coverage is limited to the near-IR, while mid-IR ob- 
servations are hampered by their lower spatial resolu- 
tions. In the optical, only the AXPs 4U 0142+61 and 
IE 1048.1-5937, and now SGR 0501+4516, were de- 
tected, but only the former at wavelengths bluer than the 
/ band. In all cases, the source of spectral information 
still relies on multi-band photometry. This comes with an 
important caveat since magnetars are known to be vari- 
able in the optical/nIR and thus the available multi-band 
photometry observations, which are often taken at dif- 
ferent epochs, are not directly comparable. Due to their 
faintness, no optical/nIR spectroscopy observations have 
been performed for any magnetar, so far. The only ex- 
ception is IE 1540.0—5408 for which both spectroscopy 
and polarimetry observations were recently performed to 
foUow-up the identification of the nIR counterpart (Is- 
rael et al. 2009a). However, at the time of writing of this 
manuscript, the data analysis is still under way. Like for 
the RPPs, useful information come from timing and po- 
larimetry observations. So far, pulsations at the X-ray pe- 
riod were detected in the optical only for the AXPs 4U 
0142+61 (Kern & Martin 2002; Dhillon et al. 2005) and 
IE 1048. 1 -5937 (Dhillon et al. 2009). No pulsation was 
detected in the nIR for 4U 0142+61 down to a pulsed 
fraction of 17% (Morii et al. 2009). Phase-averaged po- 
larimetry observations have been performed in the Kg 



band with the VLT only for IE 1048.1-5937 and XTE 
J1810-197 (Israel et al. 2009b, in preparation). 

Thus, due to the lack of a well-characterised SED, the 
origin of the optical/nIR emission of the magnetars is 
very much debated. While for both RPPs and XDINSs 
the UVOIR emission comes from the neutron star, for the 
magnetars it is still not clear yet whether the optical/nIR 
emission also comes from the neutron star. Interestingly, 
Mignani et al. (2007c) showed that, if powered by the 
neutron star rotation, the magnetar nlR efficiency L/r / E 
would be at least two orders of magnitudes larger than 
that of RPPs. This larger nIR output might be possibly 
related to the larger magnetar magnetic fields (see, e.g. 
Fig. 3 of Mignani et al. 2007c) rather than to the neutron 
star rotation, or to reprocessing of the X-ray radiation in 
a fallback disc around the magnetar. For the two magne- 
tars with the broadest optical/nIR spectral coverage, 4U 
0142+61 and IE 1048.1-5937, the SED can be fitted 
by a PL (e.g., Wang et al. 2006; Durant & van Kerkwijk 
2005), which would favour a magnetospheric origin of 
the optical/nIR emission. On the other hand, the detec- 
tion of the AXP 4U 0142+61 in the mid-IR (Wang et 
al. 2006), where it features a clear spectral bump with re- 
spect to the PL which fits the optical/nIR fluxes, has been 
considered a possible evidence for the presence of a disc. 
A similar evidence was recently found also from mid-R 
observations of IE 2259+585 (Kaplan et al. 2009). No 
other AXP/SGR has been detected in the mid-IR so far 
(Wang et al. 2007a). For IE 1048.1-5937 (Wang et al. 
2008), the deep Spitzer upper limits tend to rule out the 
presence of a disc, at least assuming the same size and 
geometry of the disc claimed around 4U 0142+61. 



The comparison between the pulsed fractions and the 
relative phases of the optical and X-ray light curves pro- 
vides an important diagnostic to discriminate between 
magnetospheric and disc emission. In the case of 4U 
0142H-61, Dhillon et al. (2005) found that the optical 
pulsed fraction (29% ± 8%) was ~5 times larger than 
the X-ray one and found no evidence for a significant 
optical-to-X-ray pulse lag, which would argue against 
the optical emission being due to X-ray reprocessing 
in a disc. However, the optical and X-ray observations 
were not simultaneous and the variation of the X-ray 
pulsed fraction along the source state might have ham- 
pered the comparison. On the other hand, in the case of 
IE 1048. 1 -5937 DhiUon et al. (2009) found, through si- 
multaneous optical and X-ray observations, that the op- 
tical pulsed fraction (21% ±7%) was actually a factor 
~ 0.7 lower than the X-ray one but they also found a 
possibly significant evidence of an X-ray-to-optical pulse 
lag (0.06±0.02). All in all, timing observations tend to 
favour a magnetospheric origin of the pulsed optical 
emission. However, a disc origin of the magnetar opti- 
cal/nlR emission is still considered a possibility (see, e.g. 
Ertan et al. 2007), at least to explain the unpulsed emis- 
sion component. 

A magnetospheric origin of the nIR emission is sup- 
ported by the erratic nIR-to-X-ray variability observed 
in XTE J1810-197 (Testa et al. 2008), which is hardly 
compatible with X-ray reprocessing in a disc. Possible 
evidence might also come from the nIR-to-X-ray decay 
rate of IE 1048.1-5937 in the 2.5 month following the 
March 2007 outburst. VLT data (Israel et al. 2009c, in 
preparation) show evidence for a larger decay in the nIR 
with respect to the X-rays, which also argues against disc 
reprocessing. The upper hmit of ~ 25% on the phase- 
averaged nIR polarisation of IE 1048.1—5937 (Israel et 
al. 2009b), with respect to the ^^10% measured in the op- 
tical for RPPs (see, e.g. Slowikowska et al. 2009 and ref- 
erences therein), does not provide compelling evidence 
for a non-thermal origin of the magnetar nIR emission, 
though. 

ccos 

For the CCOs, deep optical/nlR observations have been 
performed for nearly all sources using the HST, the VLT, 
and the Gemini. For the CCO in PKS 1209-51, the orig- 
inally proposed nIR identification with a low-mass M- 
star (Pavlov et al. 2004) was ruled out by an improved 
astrometry analysis (Mignani et al. 2007b; Wang et al. 
2007b) and by the upper limit on the proper motion of 
the proposed candidate counterpart (De Luca et al. 2009). 
No viable candidate counterpart was found for the CCOs 
in Cas A (Fesen et al. 2006), Kesteven 79 (Gotthelf et 
al. 2005), G347.3-0.5 (Mignani et al. 2008c), Puppis 



A (Wang et al. 2007b; Mignani et al. 2009d). For all of 
them, the derived optical/nIR upper hmits enable to rule 
out the presence of a binary companion other than a very 
low mass star (M5 type or later), although they do not ex- 
clude the presence of a debris disc as well, like in the case 
of the magnetars. Only for the CCO in Vela Jr. a possi- 
ble nIR counterpart (// ~ 21 .6;Ks ~ 21 .4) was identified 
by the VLT (Mignani et al. 2007d), whose nature is yet 
undetermined and could be compatible with either being 
the neutron star itself, a low-mass M-star, or a fallback 
disc. Interestingly, the Vela Jr CCO features an optical 
nebulosity detected in the R band with the VLT and also 
detected in the Ha line, which has been interpreted as a 
possible bow-shock, produced by the neutron star mo- 
tion in the ISM or a photo-ionisation nebula (Mignani et 
al. 2007b, d). However, new HST observations (Mignani 
et al. 2009e) apparently contradicts both scenarios. For 
the CCO in RCW 103, often suspected to be in a binary 
system because of its transient X-ray emission and 6 its 
hours periodicity, a candidate counterpart was proposed 
(Pavlov et al. 2004) with an M-star identified close to the 
Chandra position. However, a systematic re-analysis of 
the astrometry of all the available Chandra observations 
of the CCO (De Luca et al. 2008) ruled out the associ- 
ation with the proposed counterpart. A search for corre- 
lated X-ray and nIR variability, both along the source 6 
hours period and along a time base line of years, from all 
possible counterparts detected around the Chandra posi- 
tion was carried out by De Luca et al. (2008) using all the 
available HST and VLT data sets. However, no significant 
evidence of variability at any time scale was found. 

RRATS AND HBRPS 

For the RRATs, a search for bursting optical emis- 
sion from J1819-1458 were performed by Dhillon et 
al. (2006) but with negative results. More recently, 
nIR observations were performed with the VLT for all 
RRATs which are detected in X-rays, i.e. J1317— 5759, 
J1819-1458, and J1913H-1333. However, no counterpart 
was singled out based on colours and/or variability (Rea 
et al. 2009). The only HBRP observed in the nIR is 
PSR Jl 119— 6127, with the aim of finding evidence of 
fallback disc which might have exerted an extra torque 
on the pulsar, thus modifying its spin-down parameters 
and increasing the inferred magnetic field value. How- 
ever, VLT observations failed in detecting nlR emission 
from the pulsar position (Mignani et al. 2007c) and did 
not rule out the presence of a disc with accretion rate 
M < 3 X 10^^ g s"\ still capable of producing a substan- 
tial torque on the pulsar. 



TABLE 1. UVOIR identifications of INSs: RPPs (first group), XDESfSs (second), and magnetars 
(third). Each group is sorted according to the spin down age. Column 1 gives the INS name, while 
columns 2 to 5 give the spin down age (in logarithmic units), magnitude, distance, and interstellar 
extinction Ay. For RPPs and XDlNSs, magnitudes refer to the V band, unless otherwise indicated 
by the superscript, while for magnetars they refer to the Kj band. Column 6 to 8 give the broadband 
photometry coverage in the near-ultraviolet (nUV), optical (O), and near/mid-infrared (nIR,mIR) 
spectral bands, and a flag (Y) for spectroscopy, polarisation, and timing measurements. 
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mag 


d(kpc) 


Ay 


Phot. 


Spec. 


Pol. 


Tim. 


Crab 


^ in 


1 6 6 
1 u.u 




1 fi 

1 .U 
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Y 


Y 


R1 'sflQ— S8 


^ 1Q 


95 7* 






n 




Y 




B0540-69 


3.22 


22.0 


49.4 


0.6 


o 


Y 


Y 


Y 


Vela 


4.05 


23.6 


0.23 


0.2 


nUV, 0, nIR 


Y 


Y 


Y 


Geminga 


5.53 


25.5 


0.07 


0.07 


nUV, O, nlR 


Y 




Y 


B0656+14 


5.05 


25.0 


0.29 


0.09 


nUV, O, nlR 


Y 


Y 


Y 


B1055-52 


5.73 


24.9^ 


0.72 


0.22 


nUV, O 








B1929+10 


6.49 


25.6^^ 


0.33 


0.15 


nUV 








B1133+16 


6.69 


28 


0.35 


0.12 











B0950+08 


7.24 


27.1 


0.26 


0.03 


nUV,0 








J0108-1431 


8.3 


26.4^^ 


0.2 


0.03 


O 








J0437-4715 


9.20 




0.14 


0.11 




Y 






RX J1308.6+2127 


6.11 


28.6* 


<1 


0.14 


O 








RX J0720.4-3125 


6.27 


26.7 


0.30 


-0.3 


nUV.O 








RX J1856.5-3754 


6.60 


25.7 


0.14 


0.12 


nUV, 


Y 






RX J1605.3+3249 




26.8* 


<1 


0.06 


O 








RBS 1774 




27.2* 


0.34 


0.18 











RX J0420.0-5022 




27.5* 


0.35 


0.07 


O 








IE 1547.0-5408 


3.14 


18.5 


9 


17 


nlR 








SGR 1806-20 


3.14 


20.1 


15.1 


29 


nlR 








IE 1048.1-5937 


3.63 


21.3 


3.0 


6.10 


0, nlR 




Y 


Y 


XTE J1810-197 


3.75 


20.8 


4.0 


5.1 


nlR 








SGR 0501-^451 


4.10 


19.1 


5 


5 


0, nlR 








4U 0142-f61 


4.84 


20.1 


1.73 


1.62 


0, nIR,mIR 






Y 


IE 2259-1-586 


5.34 


21.7 


3.0 


5.7 


nIR,mIR 









SUMMARY AND CONCLUSIONS 

UVOIR observations not only add a tile to complete 
the picture of the multi-wavelength phenomenology of 
different classes of INSs but they also play a major role 
in studying the intrinsic properties of neutron stars and 
in shedding light on their formation and evolution. 

Like it is shown in the case of RPPs and XDINSs, op- 
tical/nUV observations are important, together with the 
X-ray ones, to build the thermal map of the neutron star 
surface and to investigate the conductivity in the neutron 
star interior. nUV observations are especially important 
to study the surface thermal radiation from neutron stars 
much older than 10^ years and to determine the temper- 
ature of the bulk of the neutron star surface, which is too 
low to be measured in X-rays where only very small, hot 
polar caps are seen (e.g., Becker 2009). This is crucial 
to constrain the tails of neutron star cooling models an to 
investigate possible re-heating processes (e.g.. Page 2009 
and Tsuruta 2009). While timing is instrumental in secur- 
ing the INS identification, the comparison of the UVOIR 
light curves with the radio, X-ray, and y-ray ones is im- 



portant to locate different emission regions in the neu- 
tron star magnetosphere. Moreover, the study of the giant 
pulses in RPPs, so far observed only in the radio and in 
the optical bands (Shearer et al. 2003), allows to uniquely 
study the relation between coherent and incoherent emis- 
sion. Optical polarisation measurements represent im- 
portant tests for neutron star magnetosphere models, they 
allow to measure the neutron star magnetic field and 
spin axis angles, and to unveil magneto-dynamical in- 
teractions between the neutron star and the surround- 
ing PWNe (see Mignani et al. 2007e). IR observations 
proved crucial to search for fallback discs around neu- 
tron stars, to verify accretion scenarios, to test magne- 
tars and ceo models, and to investigate the neutron star 
evolution in the immediate post-supernova phases. Fi- 
nally, UVOIR astrometry is still the best way to measure 
proper motion and parallaxes of radio-silent INSs. More- 
over, high-resolution UVOIR imaging allows to study the 
morphology and evolution of PWNe, to search for bow- 
shocks produced by the INS motion through the ISM 
and, thus, to reconstruct its spatial velocity and trace back 
its birth place and its parental stellar population. 



Table 1 summarizes the current UVOIR identification 
status for all classes of INSs, i.e. RPPs, XDINSs, and 
magnetars. The table also summarises the complete ob- 
servational data base for all the listed INSs, including 
multi-band photometry, low-resolution spectroscopy, po- 
larimetry and timing. As it is seen, a total of 25 INSs 
have been identified so far, which represents a factor of 3 
increase with respect to the number of identifications ob- 
tained by the end of the 1990s. However, this is still far 
from the number of INSs identified in the X-rays (e.g., 
Becker 2009) and from the number of INSs which are 
being identified in 7-rays by Fermi (Abdo et al. 2009, in 
preparation). Moreover, spectroscopy, polarimetry, and 
timing observations stiU represent a severe challenge for 
the current generation of 8m-class telescopes. This chal- 
lenge can be faced with the new generation of large tele- 
scopes, Uke the 42m european Extremely Large Tele- 
scope (E-ELT). Observations with the E-ELT will open 
a new era in the UVOIR astronomy of INSs, yielding to 
at least a hundred of new identification and allowing to 
carry out studies so far limited to a handful of objects 
only. 
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